Dynamics of Dzyaloshinskii domain walls in ultrathin magnetic films 
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We explore a new type of domain wall structure in ultrathin films with perpendicular anisotropy, 
that is influenced by the Dzyaloshinskii-Moriya interaction due to the adjacent layers. This study is 
performed by numerical and analytical micromagnetics. We show that these walls can behave like 
Neel walls with very high stability, moving in stationary conditions at large velocities under large 
fields. We discuss the relevance of such walls, that we propose to call Dzyaloshinskii domain walls, 
for current-driven domain wall motion under the spin Hall effect. 



I. INTRODUCTION 

Ultrathin magnetic films where the magnetic easy axis 
reorients to the film normal are one of the first systems 
in which the specificities of Nanomagnetism became visi- 
ble. A surface contribution to the magnetic anisotropy of 
a sample was predicted early by L. Neeli, and much later 
observed to dominate over the volume terms (anisotropy 
and magnetostatic) in films of thickness below 2 nm, typ- 
ically, when such samples could be fabricated. 

As the magnetic degrees of freedom along the film nor- 
mal correspond to extremely high energy modes, such 
samples were assimilated to the practical realization of 
the model system in which the simplest magnetic domain 
wall (DW) exists, namely the Bloch wall 2 -! 3 -. In this wall, 
the magnetization rotates as a spiral when one travels 
perpendicularly (direction x) across the DW, avoiding 
magnetostatic charges as divm = dm x /dx = 0, where 
rfi is the unit vector defining the magnetization direc- 
tion. The DW profile is thus governed by the exchange 
and anisotropy energies, the latter being in the ultrathin 
limit an effective anisotropy incorporating the perpen- 
dicular demagnetizing energy. As the typical DW widths 
are of a few nanometers, observations that the DW in 
such samples are of this type are extremely rare: spin- 
polarized scanning tunneling microscopy would require 
to control the orientation of the magnetic moment of the 
last atom, for example. The most convincing proof was 
recently obtained in nanostrips of (Co/Ni) multilayer by 
anisotropic magnetoresistance measurements, where for 
small nanostrip widths the DW structure was observed^ 
to transform into a Neel wall (NW), as expected from 
micromagnetics^. In the NW, the magnetization rotates 
in a cycloidal mode, i.e. with the DW normal x contained 
in the magnetization rotation plane. For perpendicular 
magnetization films that are not patterned into nanos- 
trips, the BW does not have a magnetostatic energy but 
the NW does, so that indeed the BW is the lowest energy 
DW structure, justifying the common belief. 

However, it was soon recognized that the energy den- 



sity difference between NW and BW decreases as the 
film becomes thinner—, because the 'demagnetizing fac- 
tor' of the NW reduces as the ratio iln2/(7rA) of film 
thickness t to DW width ttA. A practical consequence 
of this is the low Walker field and maximum velocity of 
the BW in these films 6 . Indeed, the stability of a DW 
structure against precession around a field applied along 
the magnetization of the domains is what limits the sta- 
tionary regime of DW motion, where the largest DW mo- 
bility (velocity to field ratio) and velocity are obtained. 
Another direct proof of this low stability of the BW in 
ultrathin films was recently obtained through the obser- 
vation, by ballistic electron emission microscopy 7 -, that 
in coupled films superimposed walls adopt an antiparal- 
lel NW structure. The stray field from the domains in 
one layer, despite its smallness because of the ultrathin 
film thickness, indeed overcomes the NW demagnetizing 
field in the other layer. 

All these considerations rest on energy terms that are 
usually considered in micromagnetics. Recently, another 
term has been experimentally uncovered in the mag- 
netism of mono- and bilayers, namely the Dzyaloshinskii- 
Moriya interaction (DMI) 8 ~ 12 . This antisymmetric ex- 
change interaction was introduced for low-symmetry 
crystal s 13 ' 14 . It was shown to favor non-uniform mag- 
netic structures 1 ^, prominently the nowadays so-called 
skyrmions. In the context of ultrathin films, A. Fert first 
mentioned that, the symmetry being reduced at the inter- 
face, such a term was also allowed^ 6 .. This was confirmed 
by an extensive analytical study 17 that considered the 2- 
sitei^ and 3-site 18 mechanisms. Quantitative calculations 
of this interaction by ab initio techniques^ confirmed its 
importance for ultrathin films. Starting with the case 
of 2 monolayers Fe epitaxially grown on the W(110) sur- 
face of a single crystal, studied in depth by spin-polarized 
scanning tunneling microscopy 20 , it was shown that the 
DMI changes the nature of the magnetic DW, forcing 
the existence of Neel walls^. The DW structure has also 
been investigated for general values of the parameters 
and orientation of the Dzyaloshinskii vector—, but in a 
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simplified model that does not solve the magnetostatic 
problem inherent to a Neel wall. Moreover, the practical 
consequences of the introduction of the DMI on the DW 
dynamics have not been investigated so far, with the ex- 
ception of one work 2 - for a soft magnetic material with 
a DMI of bulk symmetry. 

In this paper, we consider an ultrathin film with per- 
pendicular easy axis, grown on a substrate with a capping 
from a different material so that the structural inversion 
symmetry is broken along the film normal (z axis), the 
typical sample being Pt/Co(0.6 nm)/A10x^. The film 
is patterned into a strip, elongated along the x axis, of 
a width below the micrometer but not as small as to in- 
duce by magnetostatics the transition to a Neel wall^^, 
providing a well-defined geometry to study the DW prop- 
agation. As the film is ultrathin, only the x-y degrees of 
freedom exist. The conventional micromagnetic energy 
density of the system is 
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In this expression, A is the (isotropic) exchange constant, 
K u the uniaxial anisotropy constant with z the easy axis 
(sum of volume and surface anisotropics), M s the spon- 
taneous magnetization, H a is the applied field and Hd 
the demagnetizing field computed from the magnetiza- 
tion distribution through the magnetostatics equations 24 . 
In this geometry, we consider an interface DMI that reads 
in continuous form 
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This form corresponds to a sample isotropic in the plane, 
where the Dzyaloshinskii vector for any in-plane direction 
u is Dz x u with a uniform constant D, originating from 
the symmetry breaking at the z surfac o 16 ' 25 . This inter- 
face DMI differs from the bulk-like form where D[u] / / u. 
The expression ([2]) shows that indeed the DMI favors 
NWs of a given chirality, fixed by the sign of D. 

The previous studies of the DW structure under inter- 
face DM I 10 i 21 have only considered the statics (H a = 0) 
and neglected the complex effects of the demagnetizing 
field: either Hd — or a global demagnetizing effect 
was incorporated into an effective anisotropy. Here, we 
lift these two restrictions. Micromagnetic calculations 
are performed either in 2D (the full problem), or in a 
model reduced to ID (along the x axis) but with a cal- 
culation of the demagnetizing field. We also show how 
part of the physics of this DW can be captured using the 
well-known collective coordinates approach (the so-called 
q — $ model), with suitable modifications. Although most 
of the paper deals with field-driven motion of DW under 
DMI in ultrathin films, we finally discuss the relevance 
of the DMI to current-driven DW motion. 



The micromagnetic parameters chosen for the simula- 
tions are those of the typical sample considered, namely 
M s = 1100 kA/m, A = 16 pJ/m, K u = 1.27 MJ/m 3 . 
This leads to a DW width parameter A = yj A/Kq — 
5.65 nm where K = K u — ^i M^/2 is the effective 
anisotropy including the perpendicular demagnetizing 
field effect in the local approximation. For the magne- 
tization dynamics, we use the gyromagnetic ratio of the 
free electron 70 = 2.21 x 10 5 m/(A.s), and the damping 
factor a — 0.5 that is the typical experimental valued. 



II. ONE-DIMENSIONAL MODEL 

In a ID model, the degrees of freedom transverse to 
the strip are quenched. As a result, the y derivatives 
in Eqs. (P), ^ disappear. Given the typical strip- 
width (w — 500 nm) compared to the sample thickness 
(t = 0.6 nm), we decide to neglect the y component of 
the demagnetizing field. In a first approximation, justi- 
fied by the fact the nanostrip and DW widths are much 
larger than the sample thickness, we express the nor- 
mal component of the demagnetizing field by the local 
approximation Hd z — —M s m z . The longitudinal com- 
ponent of the demagnetizing field is calculated by direct 
summation with interaction coefficients C x analytically 
calculated (note also that the z component of the demag- 
netizing field can be calculated similarly, with coefficients 
C z applied to m z ). As these interaction coefficients de- 
pend on the strip width w, this value has to be specified, 
even if its influence is very small. For the numerical re- 
sults shown here, it was w = 500 nm. The numerical 
mesh size was 1 nm, and the length (in the x direction) 
of the calculation box was 300 nm. The box was shifted 
during the dynamic calculations in order to remain cen- 
tered on the DW. 

The main characteristics of the converged DW struc- 
tures, as a function of the DMI parameter D, are shown 
in Fig.[TJ The two DW moments (the integrated in-plane 
magnetization components, in nm) are plotted in Fig.[TJi. 
They show that the transition from the BW at D = 
to the NW at D > D c ss 0.13 mJ/m 2 is progressive. As 
we shall see below, this is an effect of the competition 
between the DMI and the magnetostatic energies. The 
variation of the DW energy, plotted in Fig. [1}d, conforms 
with the analytic expectation E = Eq — ttD^ at large D. 
Extrapolating the DW energy to zero, we see that this 
sample remains in the up-down state with domain walls 
up to -Dmax ~ 3.7 mJ/m 2 , above which the nonuniform 
state develops. 

Fig. [TJ; shows the variation of the DW width param- 
eter At, using the Thiele definition of this parameter 
that directly applies to the DW dynamics 2 ^. This pa- 
rameter varies weakly, but the variation is a hint that an 
analysis with a Bloch wall profile is too simple. In or- 
der to see this, the profile of the x (Neel) magnetization 
component is plotted in Fig. QJl, and compared to the 
Bloch wall profile with this parameter A. The log scale 
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stresses that, although the Bloch wall profile is very ac- 
curate around the center of the wall, it does not describe 
the so-called tails that extend far from the DW. These 
tails, well-known in the in-plane case 2 ^, are a direct con- 
sequence of the magnetostatic field inherent to the NW 
structure. Here, they are responsible for the increase of 
the wall width parameter At as the DW transforms to 
the NW, opposite to the naive expectations based only 
on the DW core energies. 

We now study the dynamics of these DW structures 
under a field H a applied along the easy (z) axis. As 
usual, a stationary motion with a high mobility is found 
at low field, followed by a lower velocity regime with a 
region of negative mobility just above the critical field 
called Walker field. Fig. shows the v(H a ) curves com- 
puted for various values of the DMI micromagnetic pa- 
rameter D (in the precessional regime, as the changes 
are periodic, the average velocity is shown). The main 
effect of DMI is to extend the stationary regime, with 
no change of mobility. The curves are very similar to 
those of the standard q — $ model for Bloch walls 28 . In 
fact, comparing to the curves of that model with identical 
Walker field, Walker velocity and damping, some differ- 
ences can be seen only in the intermediate field region 
with marked curvature. The value of the Walker field, 
plotted in Fig. (2b, shows a quasi linear increase with D. 

The quantitative effect is very large: a value D = 
0.6 mJ/m 2 increases -ffw by ten times compared to D = 
0. The Walker velocity follows exactly the same trend, 
simply because the mobility in the stationary regime is 
constant (1.96 m/(s.mT)). Such a strong effect is physi- 
cally intuitive: the DMI pins the DW magnetization to 
the Neel orientation, suppressing the precession of DW 
magnetization and thus extending the stationary regime 
to higher fields. 

The robustness of the above results was tested by com- 
parison to exact, but more demanding, 2D calculations 
(with suitably modified homemade 2 ^ and OOMMF 30 
codes), performed for different nanostrip widths w < 
500 nm and D < 0.5 mJ/m 2 . As soon as w was larger 
than the BW-NW magnetostatic transition^, the evo- 
lution of the results with w was minor. Neither a contin- 
uous rotation of the DW magnetization along the length 
of the wall, nor a tilt of the wall, were observed, so that 
the ID model is justified, as it describes a stable struc- 
ture. All quantities (DW width, Walker field, etc.) were 
found to be slightly modified, mainly because of the as- 
sumption of a local z demagnetizing field with N z — 1, so 
that the effective anisotropy is larger and the DW width 
lower. The small variations, observed above the Walker 
field, are due to the transverse degrees of freedom (see 
Ref.— for an example). They are more apparent as w or 
D increase. In order to give an idea of the magnitude 
of the 2D effects, one such calculation result is added in 
Fig. Hi. 
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FIG. 1. Equilibrium DW structure according to the ID 
model as a function of the magnitude of the DMI parameter 
D. (a) integrated DW magnetic moments for the x (Neel) 
and y (Bloch) components (big dots show the results of 2D 
calculation); (b) DW energy; (c) DW width according to the 
Thiele definition. The magnetization profile is plotted in (d) 
for the case D = 0.2 mJ/m 2 ; it shows very little variation for 
larger values of D. Curves in (a,b) were obtained with the 
q — $ model. 



III. COLLECTIVE COORDINATES q - $ MODEL 

Despite the complexity of the NW profile (Fig.QJl), and 
of the local form of the DMI, a simplified model based 
on only two collective coordinates, namely the DW po- 
sition q and the DW magnetic moment angle $, can be 
constructed, similarly to what has been done for the con- 
ventional DW 2 — . The starting point is to express the DW 



4 



120 




(b). 



120 
100 
80 



40 60 
6 z (mT) 



100 



• 1D calc. 
— q-O model 




(C) 20 
I 15 



0) 

1 10 



0.2 0.4 0.6 

D (mJ / m 2 ) 



0.8 



£ 
o 

£ 

Q 



1 1 

- X 


1 1 1 1 1 1 1 


i 1 


y y 


• IDcalc. \ 
— q-O model 






D = 0.4 mJ/m 2 

i.i.i. 





10 20 30 40 
H H a (mT) 



50 60 



FIG. 2. Field- induced domain wall motion for different value 
of the DMI parameter D, as calculated by the continuous ID 
model, (a) Velocity-field curves. The results of the 2D calcu- 
lation for a width w = 500 nm and D — 0.4 mJ/m 2 are in- 
cluded for comparison (big dots), (b) Variation of the Walker 
field Hw with D. (c) Evolution of the two DW magnetic mo- 
ments with field, in the stationary regime for D — 0.4 mJ/m 2 . 
The curves in (b,c) show the results of the q — $ model. 



surface energy tr— . As seen above, for the typical system 
that we consider the DW width is practically constant. 
We therefore ass ume th at the wall width parameter is a 
constant A = yj A/Kq. We also make the assumption 
that the magnetization rotates in a vertical plane, ro- 
tated by an angle $ around the z axis ($ = for the 
NW). The ID calculations show that this is not the case, 
neither far from the DW as m x decays slowly whereas m y 
does not, nor around the center of the wall when driven 
by a field. Nevertheless, this assumption allows a crucial 
simplification and the final results match quite well with 
ID calculations, as shown below. In these conditions, the 
DW energy reads 



a = 2AK cos 2 $ - ttD cos $ + C s 



(3) 



where K is the magnetostatic 'shape' anisotropy that fa- 
vors the Bloch wall, related to the 'demagnetizing coef- 
ficient' N x of the wall by K = N x (i M* /2. The Slon- 
czewski equations derived from Eq. ([3]) are then, with 
H K = 2K/{^M S ) and H D = ttD/ (2fi M s A), 

$ + aq/A = 7o ff a (4) 
q/A - a<I> = 70 {—Hie sin $ cos $ + H D sin $) (5) 

Minimizing a with respect to $ at H a = solves the 
statics. The stable solution is 



cos$ =ttD/ (4AK) for tt\D\ < AAK 
= sign(D) for tt\D\ > AAK. 



(0) 



Thus a critical value D c = AAK/tt appears. When |D| > 
D c we have a NW, and below this limit the DW moment 
reorients smoothly to the Bloch orientation. Figs. [T^,,b 
show the good overall agreement of this model with the 
numerical results of the ID model (for all calculations 
with the q — $ model, we used A = 5.64 nm and N x = 
0.0224, derived from the ID calculations). 

The model is also analytically solvable in the station- 
ary regime of the dynamics. Taking $ as a parameter, 
one has indeed 



H a = a sin <f> (H D - H K cos $) 



(7) 



Stability of $ requires also that Hk cos 2$ — Ho cos $ < 
0. The DW velocity v is then given by the stationary 
relation v = (70 A/a) H, where. Fig. [5t shows how well 
this reproduces the numerical results of the ID model 
for the case D = 0.4 mJ/m 2 . The highest field applica- 
ble, the Walker field, is found to be given by cos$w = 
(5 - s/6 2 + 8) /4 where 5 = D/D c = H D /H K (the ex- 
pressions for the Walker field i/w an d the Walker veloc- 
ity % are cumbersome but straightforward to write; at 
large D one has simply Hy/ « aHr>). The comparison 
with the Walker field data from the numerical solution 
of the ID model is depicted in Fig. [2Jd, again with a very 
good agreement. In the precessional regime, the dynamic 
equation for $ alone, derived from Eqs. (H]),©, should 
be integrated in order to get the oscillation period as 
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function of field. Comparison with the curve of the stan- 
dard model (D = 0) shows some difference only in the 
regime where H a is not so much above i?W; as may be 
anticipated. 



IV. COMPARISON TO EXPERIMENTS 

A first hint that the DW in ultrathin Pt/Co 
(0.6)/AlOx films is not a simple BW is provided by the 
DW dynamics under field only in this sample^ 3 -, com- 
pared to those of a series of Pt/Co(t)/Pt with t = 
0.5 — 0.8 nm 26 : the velocity in the flow regime is in- 
deed at least 5 times larger in the former case (the com- 
parison is not so direct as the Co/AlOx interface gives 
more interface anisotropy than the Co/Pt interface). The 
latter study concluded that the flow regime was preces- 
sional, with a 0.3, the stationary regime being hidden 
by the DW pinning. As a result, the high mobility ob- 
served for Pt/Co/AlOx can only be obtained in the sta- 
tionary regime, but in that case the Walker field would 
be smaller than the fields applied (up to 120 mT). As 
shown by Fig. [5] and by the analytic model, this paradox 
can be lifted with D«l mJ/m 2 or larger. This value is 
reasonable: we estimate for Fe/W(110)i£ 4 mJ/m 2 , for 
Mn/W(001)2 3 mJ/m 2 , and for Fe/Ir(lll)i 2 . 1.5 mJ/m 2 , 
supposing a dilution of the DMI values found in mono- 
or bilayers over a 3 monolayers film. 

One test of the NW structure due to DMI is to mea- 
sure the DW dynamics under a z field, with an additional 
in-plane field along the x direction. Indeed, in simple 
terms, this field favours opposite chiralities for up-down 
and down-up walls, competing with the DMI. The results 
of the full ID model for this case are shown in Fig.[3^, for 
an up-down wall with D = +0.8 mJ/m 2 . The positive 
x fields add to the DMI for stabilizing the DW moment 
in the x orientation, and increase the Walker field (and 
Walker velocity). In the expression of the DW energy c, 
the x field would enter like a DMI term, but the magneti- 
zation rotation that it induces in the domains should also 
be taken into account. This stresses the main interest of 
the DMI: it provides a large x field that exists only within 
the wall. The figure thus shows that D = 0.8 mJ/m 2 is 
equivalent to a stabilizing field (1qH x rs 128 mT, which is 
exactly the value D/(M S A) of the average over the DW 
of the equivalent B x field deduced from Eq. ©. This 
is very different from the case of a transverse (y) field, 
where small changes of mobility and Walker field are ob- 
served, depending on the sign of the z field (not shown). 

We proceed now to current-induced DW motion in 
ultrathin films with perpendicular anisotropy. Recent 
wor k 33 ' 34 , has revealed that the nature of the DWs in 
such samples is very important. Indeed, to explain 
the anomalously large DW velocity under current in 
Pt/Co/AlOx, a magnetic effect of the Rashba field at the 
interfaces has been proposed 3 ^. On the other hand, a spin 
Hall effect (SHE) due to the current flowing in the Pt or 
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FIG. 3. Dynamics of a Dzyaloshinskii DW as calculated by 
the ID model: (a) field-induced motion under an additional 
in-plane field along the nanostrip, for D — 0.8 mJ/m 2 ; (b) 
current-induced motion due to the spin Hall effect in the un- 
derlayer, for a Hall angle 9h = 0.1 rad and several values of 
D, with results of 2D calculations superposed for one case (ar- 
row depicts D increase, and curves show the large D solution 
of the extended q — $ model). 



Ta buffer layer has also been considered 3 ^. The analysis 
of the action of these torques^ 3 - shows that for perpendic- 
ular materials the SHE can efficiently move NWs at zero 
field and in opposite directions for opposite chiralities. 
This has been directly proven in the Pt/Co/Pt case 3 ^, 
by preparing NW with moderate x fields and observing 
a current assisted depinning only for the NW. In that 
structure, the sign of the SHE was moreover controlled 
by the relative thicknesses of the two Pt layers, and a 
strong DMI is not expected as the effects of both inter- 
faces compensate. We thus predict that, for asymmetric 
structures where the two interface DMI do not compen- 
sate, the DW may be a NW that efficiently moves under 
current by SHE, without the need of an applied field. 
Its direction of motion with respect to the current direc- 
tion is given by the sign of the DMI parameter and the 
SHE film position (above or below the magnetic film) . In 
addition, as the NWs induced by DMI have a fixed chi- 
rality, all walls will be pushed in the same direction under 
current by the SHE. This feature is highly desirable in 
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applications where one wants to obtain the motion of a 
train of DW in the same direction (race-track memory 
for example). It cannot be obtained with the opposite 
successive chiralities induced by an x field 3 ^. 

Calculations with the ID model (Fig.[3jD) indicate that 
DW velocity under SHE increases with D, becoming com- 
parable to experiments 23 for D > 1 mJ/m 2 for a spin 
Hall angle of 0.1 rad. The q — $ model can also be 
modified to include this term. Using the Thiele force 
equation approach, we find that the SHE can be in- 
put by an equivalent field H a = (ir/2) xM s cos $, with 
X = h9nJ/ (2e/zoM 2 i) representing the SHE-induced 
STT— , J being the current density. In the steady-state 

regime, and for S >• 1 we obtain v = vo/\Jl + (Jd/J) 2 , 
with vd = 7oAi?D and Jd = 2ateD / (ft6>#A) (note 
that the v/J slope at J = is independent of D). The 
agreement with the results of the full ID model is good 
(Fig. (3}d), despite the neglect of the magnetization rota- 
tion in the domains under the SHE. 



V. CONCLUSION 

We have studied by micromagnetics a new domain 
wall structure in ultrathin films with perpendicular 



anisotropy, as a result of uncompensated Dzyaloshinskii- 
Moriya interactions from the adjacent layers, that wc 
propose to call Dzyaloshinskii domain wall (DDW). De- 
pending on the strength D of this interaction, domain 
wall magnetizations intermediate between the Bloch and 
Neel orientation are obtained, turning to Neel at large D. 
The motion of a DDW under field is strongly affected by 
D, with a nearly linear increase of Walker field and ve- 
locity with D, up to very large values before the domains 
become unstable. Specific features of the dynamics of 
the DDW under in-plane fields have also been described. 
DDWs offer an attractive scenario to contribute to the 
understanding of the wealth of experimental results on 
currrent-induced DW motion in ultrathin films with per- 
pendicular anisotropy, and for applications. 
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